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Resumen 


La concentración de cloro como desinfectante suministrado a la red de 
distribución de agua decae con el paso del tiempo debido a su interacción 
con microorganismos, metales y otras sustancias presentes en el agua. 
Los programas de computadora, como Epanet, se utilizan para simular el 
comportamiento hidráulico y la calidad del agua en las redes de 
distribución; con su ayuda se pueden representar modificaciones físicas y 
operativas en un modelo numérico para conocer el nuevo comportamiento 
sin comprometer la calidad del servicio. La calibración del modelo de cloro 
como desinfectante consiste en tomar como referencia la concentración 
en varios puntos de la red y ajustar los coeficientes de decaimiento en el 
modelo hasta lograr que las concentraciones simuladas coincidan con las 
medidas. En el presente artículo se muestra la automatización de este 


proceso haciendo uso de la técnica heurística algoritmos genéticos y la 
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Toolkit de EPANET con MATLAB (2016). El proceso de calibración se 
realiza en un sector de la red de distribución de la ciudad de Guanajuato, 
obteniendo una correlación de 0.816 con un error promedio absoluto de 
0.08 mg/l entre las concentraciones modeladas y medidas; y que como 
resultado permitirá conocer la concentración del cloro en diversos puntos 
de la red, asegurándose de brindar un servicio de calidad durante el 
abastecimiento. 


Palabras clave: cloro libre residual, decaimiento del cloro, EPANET, 
MATLAB, Toolkit de Epanet. 


Abstract 


The concentration of chlorine as a disinfectant supplied into the water 
distribution network decays over time due to its interaction with 
microorganisms, metals, and other substances present in the water. 
Computer programs, like EPANET, are used to simulate hydraulic behavior 
and water quality in distribution networks; with them, physical and 
operational modifications can be represented in a numerical model to 
understand the new behavior without compromising the quality of service. 
The calibration of a water quality model involves referencing the 
concentration of a substance at various points in the network and 
adjusting the decay coefficients in the model until the simulated 
concentrations match the measurements. In this paper, it is shown the 
automatization of this process, using the heuristic technic genetic 
algorithms and the EPANET Toolkit with MATLAB (2016). The calibration 
process is carried out in a sector of the water distribution network of the 


city of Guanajuato, obtaining a correlation of 0.816 with an absolute 
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average error of 0.08 mg/l in the modeled concentrations with respect to 
the measurements and this will allow to know the concentration of 
chlorine at diverse points of the network, thus ensuring that quality 


service is provided throughout the supply. 


Keywords: Free residual chlorine, chlorine decay, EPANET, MATLAB, 
Toolkit. 
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Introducción 


El suministro de agua potable se realiza a través de redes que llevan el 
agua desde las fuentes de abastecimiento hasta las tomas domiciliarias. 
Este servicio debe ser continuo y con buena calidad del agua, lo que 
requiere del uso de desinfectantes para eliminar materia orgánica y 
microorganismos presentes en el agua. Uno de los desinfectantes más 
usados es el cloro, y su concentración decae con el paso del tiempo debido 
a su interacción con microorganismos, metales y otras sustancias 
presentes en el agua (Tzatchkov, Alcocer-Yamanaka, € Cortés, 2004). En 
México, la concentración de cloro en las redes de abastecimiento debe 
mantenerse dentro del rango de 0.2 a 1.5 mg/l de acuerdo con la NOM- 
127-SSA1-1994 (1994). Los programas de modelación hidráulica y de 


calidad del agua permiten a los organismos operadores simular las 
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concentraciones de cloro en las redes de distribución. El Sistema Municipal 
de Agua Potable y Alcantarillado de Guanajuato (SIMAPAG) es el 
encargado de suministrar agua potable a la ciudad de Guanajuato, 
México. Con el objetivo de ayudar a la toma de decisiones en aspectos de 
calidad de agua, se modeló la concentración de desinfectante en un sector 


de su red de distribución y con ello validar el algoritmo planteado. 


Diversos investigadores han trabajado en el modelado de la calidad 
del agua en redes de abastecimiento de agua, contribuyendo a la 
comprensión del decaimiento del cloro. Se han desarrollado diferentes 
métodos para modelar el decaimiento de una sustancia mediante el uso 
de coeficientes de reacción en el seno de fluido (kp) y en la pared de la 
tubería (Kw) (Conagua, 2016). Vasconcelos, Rossman, Grayman, Boulos 
y Clark (1997) analizaron los principales factores que influyen en el 
decaimiento del cloro en acueductos, plantas de tratamiento y sistemas 
de distribución de agua potable en EUA. Alcocer-Yamanaka, Tzatchkov y 
Cortés (2004) analizaron el decaimiento del cloro en una red de 
abastecimiento en Culiacán, Sinaloa, México. Munavalli, Kumar y Kulkarni 
(2009) modelaron la calidad del agua en una red de distribución de la 
India en estado permanente. Díaz-Duque (2015) modeló la calidad del 
agua en la ciudad de Bogotá, Colombia, usando algoritmos genéticos y 
diferentes órdenes de reacción para representar el comportamiento de kb, 
asignando un mismo valor de ky a tuberías con características similares 
(material y edad de las tuberías, dirección del flujo, etc.). Baño-Caballero 
(2016) tomó en cuenta el impacto que tiene la temperatura del agua para 
representar los valores de kp» en la modelación de calidad de agua en el 
municipio de Alcantarilla, España. Minaee et al. (2019) retomaron el 


trabajo de Munavalli et al. (2009), utilizando su caso de estudio, 
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mejorando sus resultados al implementar otras técnicas a la utilizada por 


el autor inicial. 


El programa de Epanet simula la calidad del agua con aproximación 
Lagrangiana a intervalos fijos de tiempo en segmentos discretos de 
tubería y modela el proceso convectivo sin realizar intercambio de masa 
(Rossman, 2001). Al avanzar el flujo, las sustancias disueltas pueden ir 
generando reacción a una velocidad creciente o decreciente, por lo que 
requiere de intervalos de paso mucho menores que en la simulación 
hidráulica (Hernández-Cervantes, 2021). El cloro en la red reacciona y 
decae, el decaimiento está representado por la reacción en el seno del 
flujo y por la reacción por el contacto con la pared de la tubería. La 
primera reacción se mide a través del coeficiente de reacción en el seno 
del flujo, k», y se obtiene midiendo el cloro residual libre en distintos 
intervalos de tiempo en muestras de agua de la red (Alcocer-Yamanaka 
et al., 2004), y describe la representación del grado exponencial de 


decaimiento mediante la Ecuación (1): 
0 =E20* (1) 


Donde C(t) es la concentración medida en el tiempo t, en mg/l; Co, 
la concentración inicial en mg/l; t, el instante de tiempo de medición en 


horas, y k» es el coeficiente de reacción en el seno del flujo en h”!. 


La determinación de ky es más compleja que la de k, debido a que 
Kw depende del grado de transferencia de masa de los materiales alojados 
en las paredes de la tubería y de las características de ésta, especialmente 


del material y la antigúedad. En las tuberías metálicas, la rugosidad 
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aumenta con su antigúedad, y con ello las incrustaciones y los efectos de 
los agentes corrosivos que impactan en este coeficiente de reacción. El 


coeficiente Kw está implícito en la Ecuación (2): 


KywKkf 


K =Kp + Ru(Ky+K f) 


(2) 


Donde ky es el coeficiente de reacción en pared de tubería en m/h; 
Kf, el coeficiente de transferencia de masa desde el interior de fluido hacia 
las paredes de la tubería en m/h; Ra, radio hidráulico de la tubería en m, 


y K es el coeficiente total de decaimiento, h"!, 


La representación correcta de un modelo se comprueba mediante el 
coeficiente de correlación, el cual representa la relación entre dos 
variables: valor de 1 o cercano significa que una variable aumenta o 
disminuye si la otra también lo hace; un valor de -1 o cercano significa 
que una variable aumenta si la otra disminuye, y viceversa; un valor de 
O o cercano implica que si una de las variables aumenta o disminuye, la 
otra no se ve afectada. En la modelación se busca una correlación alta, 
aunque no garantiza que el modelo esté correctamente calibrado, pues 
dos variables pueden estar relacionadas sin ser iguales. Otro parámetro 
utilizado es el error absoluto promedio de las concentraciones (Ecuación 
(3)). El objetivo es tener errores bajos, lo cual se traduce en 


concentraciones simuladas y muestreadas muy parecidas entre sí: 


NN 2 |Cm=Cs| (3) 
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Donde Eprom es el error absoluto promedio en mg/l; Cm, la 
concentración de cloro muestreada en mg/l; Cs, la concentración de cloro 


simulada en mg/l, y Med es el número de mediciones realizadas. 


Los Algoritmos Genéticos (AG) son métodos adaptativos utilizados 
para resolver problemas de búsquedas y optimización (Holland, 1975; 
Beasley, Bull, 8 Martin, 1993). Son útiles para tratar con problemas no 
lineales, parámetros no continuos y objetivos de optimización compleja. 
Fueron propuestos por Holland (1975) como un algoritmo basado en un 
tipo darwiniano de estrategia de supervivencia con reproducción sexual, 
donde los individuos fuertes de una población tienen una probabilidad alta 
de crear descendencia (Díaz-Arévalo, Izquierdo-Sebastián, López- 
Jiménez, 8 Pérez-Garcia, 2004). Se parte de un conjunto inicial de 
individuos, llamado población, generado de manera aleatoria; cada 
individuo representa una posible solución al problema (Gestal, Rivero, 
Rabuñal, Dorado, € Pazos, 2010) y tendrá asociado un valor, llamado 
factor de aptitud, que califica su validez como solución al problema. En 
función de este valor se le darán más o menos oportunidades de 
reproducción. Además, con cierta probabilidad se realizarán mutaciones 
(Goldberg, 2002). El proceso de reproducción se repite durante varias 
generaciones hasta llegar a individuos que representen una solución 
aceptable. El tamaño del problema, complejidad y número de variables 
determina el número de generaciones necesarias (Díaz-Arévalo et al., 
2004). Se implementará esta técnica para encontrar la solución óptima 
para los valores de ky que representen de mejor manera la evolución del 
desinfectante a lo largo de la red. Esta técnica ha sido utilizada 
anteriormente para la optimización hidráulica y de calidad en redes de 
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abastecimiento de agua (Ozdemir €: Erkan-Ucaner, 2005; Nouiri 8 Lebdi, 
2006; Hernández-Cervantes, 2015; Díaz-Duque, 2015; Hernández- 
Cervantes, Mora-Rodríguez, Delgado-Galván, Ortiz-Medel, 8 Jiménez- 
Magaña, 2016; Nouiri, 2017; Minaee et al., 2019; Pineda-Sandoval, 
2020), mostrando con éxito su capacidad de llegar a buenos resultados 
en problemas complejos de optimización. La presente investigación 
muestra un proceso de calibración, considerando agrupamiento de 
tuberías de diámetro y rugosidad para la obtención de Ky óptimo, con una 
función de aptitud novedosa que requiere de un número básico de 
elementos tanto en técnica heurística y computacional para proveer 


resultados adecuados de calidad al sistema de abastecimiento de agua. 


Zona de estudio 


El Sector Filtros 8-12 (SF812) se localiza en la parte norte de la ciudad de 
Guanajuato, México; se abastece del agua proveniente de dos presas; la 
red es principalmente de tipo abierta y distribuye el agua por gravedad 
(Figura 1). Comienza en la planta potabilizadora llamada *Filtros”. Al final 
de ésta se tiene un tanque de almacenamiento que suministra cinco 


sectores. 
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Leyenda 


* Tanque 


— Tuberías 
Ea Sector Filtros 8-12 


Figura 1. Red de distribución del SF812. 


La red de distribución del SF812 tiene más de 30 kilómetros de 
longitud. Con base en los diámetros y rugosidades de las tuberías, la red 


se conforma de la siguiente manera: 


e Red primaria: compuesta por tuberías de fierro fundido (FoFo). Sus 


diámetros van desde los 150 hasta los 600 mm. 
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e Red secundaria: compuesta principalmente por tuberías de PVC, con 


algunas tuberías de FoFo. Sus diámetros van desde los 13 hasta los 
125 mm. 


El consumo promedio diario en el sector es de 33.88 l/s, con un 
consumo máximo de 44.31 l/s a las 13:00 h y un consumo mínimo de 
19.3 l/s a las 4:00 horas. El modelo hidráulico está compuesto por 1 226 
nodos, un reservorio, un tanque de almacenamiento, 1 070 tuberías y 
193 válvulas (Figura 2); la calibración hidráulica que se utilizó para la 
presente investigación fue de septiembre de 2019, con un coeficiente de 


correlación de las presiones de 0.988. 
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Figura 2. Modelo hidráulico en EPANET del SF812. 
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Modelación de la calidad del agua 


Metodología 


Para modelar la calidad del agua en la red del SF812, además de tener el 


modelo hidráulico calibrado, se planteó la siguiente metodología: 


a) Medir la concentración de cloro inyectada desde la potabilizadora. 
b) Obtener el coeficiente de reacción en el seno del agua (k»). 

Cc) Medir la concentración de cloro en puntos clave de la red. 

d) Calcular la concentración de cloro en Epanet. 

e) Automatizar la calibración mediante código propuesto con AG. 

f) Verificar parámetros de AG en red teórica. 

g) Calibrar coeficientes ky en red SF812. 


h) Realizar un análisis de sensibilidad del modelo. 


a) Medir la concentración de cloro inyectada 


La concentración de cloro inyectada a la red de distribución sirve como 
valor de frontera para las concentraciones que se tendrán a lo largo de la 
red y se obtuvo mediante un controlador de cloro (Figura 3) instalado en 
la tubería principal a cinco metros de la salida del tanque de 


almacenamiento. 
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Figura 3. Controlador instalado en filtros. 


Se realizó un registro a cada hora, durante 51 días, desde el 25 de 
noviembre de 2019 hasta el 16 de enero de 2020. Con los datos 
registrados se generó un patrón de concentración diaria, el cual se 
muestra en la Figura 4. La dosis de cloro que se suministra en la última 
etapa de potabilización es constante, pero como se puede observar en el 
registro de la concentración de cloro en el punto de medición en la tubería 
principal de la red, la concentración varía a lo largo de las 24 horas. Esto 
puede deberse a las características del agua superficial y la variación de 


la temperatura a lo largo del día (Geng, Fan, Shi, Zhang, 8: Li, 2022; Hua 
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et al., 2021; Zhong et al., 2021). Además, el tanque, al ser de 


mampostería y con la turbulencia del agua en su interior, genera erosión 


en sus paredes (Angeloudis, Stoesser, € Falconer, 2014), lo cual también 


podría incrementar la reacción del cloro. 


00 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 19 20 21 22 23 
Hora 


Figura 4. Concentración de cloro que ingresa al SF812 media 
horaria del 25 de noviembre de 20109 al 16 de enero del 2020. 


b) Obtener el coeficiente de reacción k; 


Para determinar el valor del coeficiente de reacción kp se realizaron 
pruebas de frasco en dos días diferentes: el primero el 20 de septiembre 
de 2019 y el segundo el 24 enero de 2020. Ambas muestras de agua se 
tomaron a las siete de la mañana, en el tanque de almacenamiento, 


posterior a la potabilización, a una profundidad media para evitar 
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contaminación de la superficie o del fondo. Las muestras se toman en un 
frasco de vidrio y se protegieron de la luz solar, para evitar que influyera 
en el decaimiento del cloro. Las concentraciones se midieron mediante un 
fotómetro HI 96711 de la marca HANNA, el cual mide el contenido de 
cloro libre y total (Cl2) en muestras de agua con un rango de O a 5 mg/l 
y precisión de +0.03 mg/l. La Figura 5 muestra los resultados de ambas 
pruebas de frasco con el que se generó la curva de decaimiento del cloro 
debido a la reacción con las características del agua de Filtros (es el 
nombre de la potabilizadora, por eso va en mayúscula). Se aproximó la 
curva de regresión exponencial de cada prueba y se obtuvo el coeficiente 
kp correspondiente al tiempo (7) en una ecuación de primer orden. Los 
resultados de las pruebas se muestran en la Tabla 1 (Pineda-Sandoval, 
2020). 


1.6 
1.4 A e Primera prueba (septiembre de 2019) 
OS e Segunda prueba (enero de 2020) 
1.2 k | 
es. Exponencial (primera prueba, septiembre de 2019) 
1 E E A Exponencial (segunda prueba, enero de 2020) 
“0 
bd . 
.. 0 
0.6 - 0 
50 
0.4 (0, e 
Q.. [OS 0.3036-9-104"T A e, - 

o o ¡mamma lo ¿PON 0 as ¡C = 1.4238e:0.063 *T 
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A PON O O tt... rro O 
0 40009000090..O O A 2) 
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Figura 5. Decaimiento del cloro por reacción en el seno del agua. 


Tabla 1. Coeficientes k, a partir de pruebas de frasco. 


Coeficiente de reacción kp 
Prueba de frasco 
1/hora 1/día 
Primera prueba -0,063 -1.512 
Segunda prueba -0.104 -2.496 


Aun cuando las muestras se tomaron a la misma hora, pero en 
fechas diferentes, la temperatura ambiente mínima en el mes de 
septiembre oscila alrededor de los 14 *C y en el mes de enero oscila 
alrededor de los 7 “C. Además de ese factor, el sistema de agua registró 
una disminución de la concentración inyectada por cuestiones operativas. 
Ambos factores hicieron que la concentración inicial de la segunda prueba 
fuera baja respecto a la primera y, con ello, el coeficiente de reacción kb» 


obtuvo un valor diferente. 


Posterior a la toma de muestras, los frascos se llevaron al interior 
del laboratorio y en las horas subsecuentes se midió la temperatura. Para 
la primera prueba, la temperatura osciló entre 17.0 y 27.7 *C; mientras 
que para la segunda varió entre 13.0 y 24.0 *C. Al ser menor la 
temperatura de la segunda prueba, el cloro debería decaer a menor 
velocidad, como mencionan García-Ávila et al. (2020), y Hua, Vasyukova 
y Uhl (2015). Sin embargo, el coeficiente de reacción de la segunda 


prueba fue mayor posiblemente por su concentración inicial de 0.53 mg/l 
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que, como lo mencionan Fisher, Kastl y Sathasivan (2012), con 


concentraciones alrededor a 0.5 mg/l tienen decaimiento acelerado. 


Cc) Medir la concentración de cloro en la red 


Se realizaron mediciones locales del cloro libre residual en diversos sitios 
de la red. Los puntos de muestreo se seleccionaron con base en ciertas 
características, como zonas de alto consumo, puntos alejados de la fuente 
y zonas con tiempos de residencia del agua altos; estos tiempos de 
residencia altos se obtuvieron en el programa de Epanet a partir del 
modelo hidráulico validado del sector en periodo extendido. En conjunto 
con personal de SIMAPAG se realizaron campañas de medición durante 
los meses de septiembre y octubre de 2019, y febrero de 2020, 
obteniendo un total de 35 sitios de medición, que representan el 2.85 % 
de los nodos que conforman al modelo, superior al 2 % recomendado para 
calibración de modelos de calidad del agua en redes de abastecimiento 
(AWWA, 2017). La Figura 6 muestra la ubicación de los sitios 


muestreados. 
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Figura 6. Localización de los sitios de muestreo en el SF812. 


d) Calcular la concentración de cloro en Epanet 
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Para calcular la concentración de cloro en Epanet es necesario llevar a 
cabo la simulación en periodo extendido, por lo que se debe fijar un 
tiempo total de simulación en el navegador de Epanet y establecer la 
calidad inicial en las fuentes; para el caso del sector SF812, se determinó 


un patrón de concentración con los datos de la Figura 4. 


Posteriormente, hay que fijar el valor del coeficiente de reacción de 
masa ko», para lo cual se utilizó el valor de la primera prueba de la Tabla 
1. Al ser un valor constante para todas las tuberías, se puede ingresar 
desde el navegador de Epanet. En opciones de reacción se inserta el valor 
global del coeficiente de reacción de masa. Además, en esas mismas 
opciones se verifica que el orden de reacción tanto de masa como de 
pared sea de primer orden, considerado adecuado para simular el 


decaimiento del cloro. 


Por último, en las opciones de tiempo del mismo navegador de 
Epanet se verifica que el paso del tiempo de calidad sea de cinco minutos, 
valor que se considera estándar para este tipo de simulaciones. Al realizar 
la simulación se pueden observar las concentraciones resultantes en los 
nodos en cada instante de tiempo. El coeficiente de reacción de masa es 
el valor a calibrar con la presente propuesta de investigación, por lo que 
en el siguiente paso se muestra cómo se genera dicho parámetro de 


manera automatizada. 


e) Automatizar la calibración mediante código 
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El programa de calibración desarrollado en el presente proyecto se 
denomina Código de Calibración de ky mediante Diámetros y Rugosidades 
(CCKwDyR). El parámetro que se utiliza para calibrar el modelo es el 
coeficiente de reacción de pared (kw), y se establece asignar el mismo 
coeficiente Kw a las tuberías que tengan las mismas características de 
diámetro y rugosidad (Díaz-Duque, 2015). El código CCKwDyR usará la 
toolkit de Epanet para obtener la información desde el modelo hidráulico 
y procesar sus resultados (Eliades, Kyriakou, Vrachimis, €: Polycarpou, 
2016). La Figura 7 muestra el diagrama de flujo del funcionamiento 
general de los AG utilizado por el CCKwDyR para obtener soluciones 


óptimas de valores de k,. 
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Figura 7. Diagrama de flujo del funcionamiento general del 
CCKwDyR. 


Con base en los valores de ky encontrados por diversos autores, se 


propone el siguiente rango de valores para Ky en el CCKwDyR: 
0.0 > kw > —1.5m/día 


Se asigna con mayor frecuencia valores de ky cercanos a 0, y con 
menor frecuencia valores cercanos a -1.50 debido al comportamiento 
exponencial del decaimiento del cloro. El rango de valores ky se presenta 


de la siguiente manera: 


e Para -0.01 > Ky > -0.50, se tendrá un incremento de 0.01. 
e Para -0.50 > Ky > -1.00, se tendrá un incremento de 0.02. 


e Para -1.00 > Ky > -1.50, se tendrá un incremento de 0.05. 


En CCKwDyR, la Función de Aptitud (FA) empleará el error por 
medición promedio para evaluar a los individuos (Ecuación (3)) y ésta 
varía en el rango de 0.5 a 1, donde los valores cercanos a 1 representan 
menor error de medición promedio de las concentraciones simuladas 
respecto a las medidas. La FA propuesta se basará en la progresión de 
una función sigmoide, donde se muestra una progresión lenta cuando el 


error por medición promedio es alto; sin embargo, cuando este error 
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comienza a ser pequeño, su progresión se vuelve rápida hacia valores de 
FA cercanos a 1. Se incluye el coeficiente de 33.22192436 con la intención 
de separar a los individuos con error grande de aquellos que presentan 
un error pequeño, considerado el punto de esa separación con un valor 
de error por medición promedio de 0.15 mg/l y un FA de 0.55 (Figura 8). 


La FA se describe en la Ecuación (4): 


FA = —,-—— (4) 


1+e 33.22192436-Eprom 


Donde FA es la función de aptitud, adimensional, y Eprom es el error 


absoluto promedio en mg/!. 
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Figura 8. Función de aptitud utilizada por el CCKwDyR. 
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f) Verificar parámetros de AG en red teórica 


Para verificar el funcionamiento del código desarrollado, así como para 
determinar los parámetros necesarios para que el CCKwDyR pueda 
encontrar una solución aceptable (número de individuos y generaciones), 
se realizó una calibración en una red teórica, llamada Net3 (Figura 9), que 
viene incluida como ejemplo en el programa Epanet, y que cuenta con los 
siguientes elementos: dos fuentes de abastecimiento, tres tanques, 97 
nodos, dos equipos de bombeo y 117 tuberías, las cuales tienen los 


siguientes diámetros y rugosidades: 


e Diámetros (mm): 203, 254, 305, 356, 406, 457, 508, 610, 762 y 2 
515, 


+ Rugosidad (Ch): 110, 130, 140, 141 y 199. 
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Figura 9. Red de distribución Net3. 


Para generar las concentraciones de cloro en la red teórica se utilizó el 
coeficiente kp promedio = -2.004 y se propusieron valores de ky dentro 
del rango de O a -0.4 m/día, con mayor decaimiento conforme disminuye 
el diámetro y la rugosidad; se registró la concentración en 17 de los nodos 
(17.5 % de los nodos totales), a manera de mediciones puntuales. El 


CCKwDyR buscó la solución óptima de ky dentro del rango de O a -1.5 
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m/día. La simulación hidráulica se realizó durante 144 horas (seis días) 
para estabilizar las concentraciones de cloro a lo largo de la red, y 
solamente se analizaron las concentraciones durante las últimas 24 horas 


de la simulación. 


Los resultados de la verificación se muestran en la Tabla 2. Se observa 
que es posible obtener resultados adecuados con un número reducido de 
individuos y generaciones, gracias a la agrupación de las tuberías con 
base en su diámetro y rugosidad, que en el caso de la NET3 se tienen 14 
agrupaciones a los que se buscará un valor de Ky óptimo, a diferencia de 
un algoritmo, donde cada tubería tiene su propio valor de ky, que para el 
caso de la NET3 sería de 117 valores de k, lo que reduce el número de 


opciones posibles. 


274 
o = 2024, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 
Open Access bajo la licencia CC BY-NC-SA 4.0 15(2), 248-303. DOI: 10.24850/j-tyca-15-02-06 


(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


Tecnología y 


CienciaszAgua 


OPEN Qucc.ss 


1) Check for updates 


Tabla 2. Evaluación de parámetros para el algoritmo genético. 


Escenario Prueba Individuos | Generaciones FA E US Correlación 
promedio (mg/l) 
1 100 10 0.9830 0.0074 0.999 
1 2 100 10 0.8881 0.0145 0.998 
3 100 10 0.9210 0.0123 0.999 
1 50 10 0.8427 0.0179 0.997 
2 2 50 10 0.9700 0.0087 0.998 
3 50 10 0.8973 0.0139 0.999 
1 100 5 0.8820 0.0150 0.999 
3 2 100 5 0.8675 0.0160 0.998 
3 100 5 0.7986 0.0219 0.996 
1 50 5 0.7484 0.0276 0.997 
4 2 50 5 0.8628 0.0164 1.000 
3 50 5 0.7708 0.0248 0.997 
1 100 20 0.9860 0.0071 1.000 
5 2 100 20 0.9801 0.0077 1.000 
3 100 20 0.9853 0.0072 1.000 
1 200 10 0.9717 0.0085 1.000 
6 2 200 10 0.9756 0.0082 0.999 
3 200 10 0.9985 0.0046 1.000 
1 200 20 0.9983 0.0047 1.000 
7 2 200 20 0.99751 0.0050 1.000 
3 200 20 0.9978 0.0049 1.000 
(cc) OSO; 275 
AMA 024, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


15(2), 248-303. DOI: 10.24850/j-tyca-15-02-06 


a 1) Check for updates 
OPEN (Y ACCESS 


Tecnología y 
CienciaszAgua 
A partir de estos resultados se determina que los parámetros para 
simular los AG en la red SF812 serán de un número de individuos de 400, 
tomando en cuenta que el número de agrupaciones de tuberías de 
diámetros y rugosidad para la red SF812 es de 30 y, con ello, el número 
de valores a buscar de ky óptimo. Por otro lado, al tener un número de 
elementos mayor que la Net3 y considerando que el tiempo de procesador 
no sea muy extenso para la calibración, se propone un número de 
generaciones de 30, verificando en la Tabla 2, que al aumentar el número 
de generaciones se mejoran los resultados. De igual forma se aumenta la 
participación de los individuos durante la aplicación del algoritmo 
genético, incrementando la probabilidad de cruza a 90 %, y se ajusta la 
probabilidad de mutación al 5 % para mantener variedad en la población. 
El alza del número de individuos y generaciones ayuda a evitar que el 
algoritmo llegue a un resultado óptimo local en lugar del óptimo global 
(Díaz-Arévalo et al., 2004); sin embargo, utilizar valores demasiado altos 
generaría tiempos de simulación excesivos, y requeriría mayor capacidad 


computacional. 


g)  Calibrar coeficientes k,, en red SF812 


Se llevaron a cabo tres pruebas para obtener la calibración del modelo, 
utilizando un tiempo de simulación de 288 horas para que el modelo 
pudiera representar el decaimiento del cloro correctamente en los puntos 
más alejados de la fuente, pues conforme se aleja del punto de inyección 
los efectos del decaimiento en las concentraciones de cloro tardan en 


hacerse presentes (Pineda-Sandoval, 2020); y solo se obtienen los 
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resultados de la calibración en las últimas 24 horas, donde las 
concentraciones de cloro en el modelo se estabilizan. En las tres pruebas 
se llegó a resultados similares, por lo que se puede asumir que el modelo 
llegó al mejor resultado posible, con un coeficiente de correlación de 
0.769 y error absoluto promedio de 0.097 mg/l (Tabla 3). 


Tabla 3. Resultados de la calibración del SF812. 


Concentraciones de cloro (mg/l) 
Parámetro Prueba 1 Prueba 2 Prueba 3 
FA 0.5779 0.5768 0.5771 
Error absoluto (mg/l) 0.097 0.098 0.098 
Correlación 0.769 0.767 0.764 


La Figura 10 muestra las concentraciones simuladas en Epanet en 
comparación con las concentraciones muestreadas en campo. Si bien se 
obtuvo un error absoluto promedio dentro del rango obtenido por otros 
autores (Tabla 1), la correlación del modelo es más baja. Este resultado 
de calibración se puede deber a diversos factores, entre los que destaca 
que las concentraciones muestreadas representan una tasa de 
decaimiento diferente entre los dos periodos de muestreo y la variación 
en los consumos de agua en el sector para las distintas fechas con las que 


se tienen de los datos. 
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O Concentraciones calibración pru 


e A Concentraciones calibración pru 
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Concentración de cloro (mg/l) 


9 


de los nodos muestreados 


Figura 10. Calibración inicial de la concentración de cloro en el 
sF812. 


h) Realizar un análisis de sensibilidad del modelo 


Se realizó un análisis de sensibilidad de los parámetros que intervienen 
en el modelado de la calidad del agua, con el objetivo de verificar su 
influencia en la simulación y encontrar posibles fuentes de error en la 


calibración del modelo. Para ello, se propuso un escenario hipotético del 


278 
EME > 024, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 
Open Access bajo la licencia CC BY-NC-SA 4.0 15(2), 248-303. DOI: 10.24850/j-tyca-15-02-06 


(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


seal Qucc.ss (Mi) Check for updates 
Tecnología y 
Ciencias¿Agua 
SF812 calibrado, considerando el patrón de inyección de cloro de la Figura 
4, Se usará el kp calculado en septiembre de 2019 y coeficientes ky con 
valores entre O y -0.4 m/día. Se determinó la concentración de cloro en 


47 nodos del modelo, simulando un muestreo en campo. 

Se evaluará la sensibilidad del modelo respecto a cinco parámetros: 
e Aumento y disminución en la inyección inicial de cloro. 
e Aumento y disminución del coeficiente ko». 


e Aumento del error en las concentraciones muestreadas utilizadas para 


la calibración. 
e Aumento y disminución del consumo de agua de la red. 


e Distribución de los puntos muestreados (omitiendo muestreo en línea 


principal, secundaria o puntos extremos). 


Se utilizaron los mismos parámetros del AG que se usaron en la 
calibración inicial (número de individuos, número de generaciones, 
probabilidad de cruza y probabilidad de mutación), evaluando así, si éstos 
son suficientes para llegar a una buena calibración en una red real como 
la del SF812. La Tabla 4 resume las variaciones propuestas para este 


análisis de sensibilidad. 
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Tabla 4. Resumen de escenarios para el análisis de sensibilidad. 


Análisis de parámetros 


Escenario Inyección | Coeficiente Error de yA 
, EN Descripción de la prueba 
(mg/l) ky (día”*) | medición (mg/l) 
Inicial 1.2318 -1.512 +0 Red calibrada inicialmente 
1 1.8477 -1.512 +0 +50 % inyección 
Inyección inicial de 2 1.3856 -1.512 +0 Mayor inyección registrada 
cloro 3 0.8334 -1.512 +0 Menor inyección registrada 
4 0.6159 -1.512 +0 -50 % inyección 
5 1.2318 -2.268 +0 +50 % kp medido 
6 1.2318 -1.890 +0 +25 % kp» medido 
Coeficiente kb» 
7 1.2318 -1,134 +0 -25 % kp medido 
8 1.2318 -0.756 +0 -50 % kp medido 
9 1.2318 -1.512 +0.03 Rango de error del controlador 
10 1.2318 -1,.512 +0.09 Error intermedio 
Errores de medición 
11 1.2318 -1.512 +0.15 Máximo error del Estado del Arte 
12 1.2318 -1.512 +0.25 Error extremo en las mediciones 
13 1.2318 -1.512 +0 +25 % consumo base 
14 1.2318 -1.512 +0 +50 % consumo base 
Consumos de agua 
15 1.2318 -1.512 +0 -50 % consumo base 
16 1.2318 -1.512 +0 -25 % consumo base 
17 1.2318 -1.512 +0 Sin mediciones en la línea principal 
18 1.2318 -1.512 +0 Sin mediciones intermedias 
Muestreo 19 1.2318 -1.512 +0 Sin mediciones en los extremos 
-50 % datos muestreados (24 
20 1.2318 -1.512 +0 


mediciones) 
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La Tabla 5 resume los resultados de correlación y error absoluto 


los resultados manteniendo 


La columna del 


análisis de 


la variación del 


parámetro, mientras que la columna del escenario calibrado presenta los 


resultados de la calibración sin variación en ninguno de los parámetros. 


Tabla 5. Resumen de resultados del análisis de parámetros. 


Análisis de sensibilidad Escenario calibrado 
Prueba FA . - 
Correlación | Error absoluto (mg/l) Correlación | Error absoluto (mg/l) 

Inicial 0.955 0.999 0.010 1.000 0.005 
1 0.721 0.965 0.032 0.966 0.172 
, 2 0.806 0.989 0.021 0.988 0.064 

Inyección 
3 0.598 0.963 0.076 0.963 0.173 
4 0.548 0.947 0.156 0.945 0.230 
5 0.885 0.996 0.015 0.990 0.108 
6 0.852 0.996 0.017 0.994 0.052 

Coef. kb 
7 0.936 0.998 0.011 0.997 0.046 
8 0.905 0.997 0.013 0.996 0.090 
9 0.712 0.988 0.032 0.994 0.024 
10 0.595 0.950 0.079 0.989 0.030 

Medición 
11 0.561 0.894 0.125 0.97 0.092 
12 0.535 0.768 0.220 0.956 0.113 
13 0.902 0.998 0.014 0.998 0.016 
14 0.873 0.995 0.016 0.998 0.022 

Consumos 
15 0.716 0.986 0.033 0.988 0.080 
16 0.875 0.997 0.015 0.998 0.020 
17 0.816 0.988 0.020 0.993 0.022 
18 0.914 0.998 0.013 0.996 0.015 

Muestreo 
19 0.936 0.998 0.011 0.998 0.011 
20 0.844 0.995 0.018 0.992 0.025 

281 


Open Access bajo la licencia CC BY-NC-SA 4.0 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


Tecnología y ciencias del agua, ISSN 2007-2422, 


15(2), 248-303. DOI: 10.24850/j-tyca-15-02-06 


a 1) Check for updates 
OPEN ACCESS 


Tecnología y 


CienciaszAgua 


Con el análisis realizado se observa lo siguiente: 


El CCKwDyR es capaz de llegar a calibraciones muy buenas cuando los 
parámetros suministrados son los correctos, lo cual valida su 
funcionamiento. Indica además que el número de individuos y 
generaciones usado (400 y 30, respectivamente) es suficiente para 


lograr buenas calibraciones en redes como la SF812. 


La concentración inicial de cloro es uno de los parámetros que más 
afecta a la calibración del modelo. Cuando se sobreestima, los nodos 
más cercanos al punto de inyección son los únicos afectados; cuando 
se subestima, las concentraciones de cloro en todo el modelo se ven 


afectadas, con valores por debajo de los muestreados. 


La variación en el coeficiente kp no es fácil de detectar y puede ser 
compensada por el CCKwDyR mediante los valores de ky propuestos 
por el algoritmo. Entre más certeza se tenga de que el valor de kp es 
el correcto, mayor probabilidad se tiene de que los valores de Kw 


representen a los valores reales y no sólo una solución matemática. 


e Un muestreo impreciso de las concentraciones de cloro afecta a la 
calibración del modelo de manera grave. Sin embargo, se observó que 
al evaluar la calibración con las concentraciones originales, el modelo 
mejoró su correlación y disminuyó el error promedio; esto sugiere que 
el CCKwDyR es capaz de compensar (hasta cierto punto) los errores 
de las mediciones. No obstante, se requieren más pruebas para 


verificar esta observación. 
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e Las variaciones pequeñas en el consumo de agua no afectan de 
manera importante la calibración (aunque se recomienda verificar la 


calibración hidráulica antes de modelar la calidad del agua). 


e Los muestreos en la línea principal son los más importantes para 
representar de forma correcta las concentraciones en toda la red. Se 
concluye además que un muestreo del 2 % de los nodos es suficiente 
para calibrar el modelo (AWWA, 2017). Sin embargo, de ser posible 


es preferible contar con un muestreo mayor. 


Resultados 


Con base en lo observado en la calibración inicial y en el análisis de 
sensibilidad, cuando se sobreestima la inyección de cloro, las 
concentraciones en los nodos más cercanos al punto de inyección no 
alcanzan a decaer lo suficiente como para obtener un buen ajuste del 
modelo. Un fenómeno similar ocurre en las calibraciones del SF812 (Tabla 
6). Los errores más grandes ocurren en las tuberías más cercanas al 
punto de inyección, mientras que en el resto de los nodos las 
concentraciones se ajustan de mejor manera, lo que sugiere que la 


inyección inicial de cloro debe ser menor para lograr un mejor ajuste. 
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Tabla 6. Concentración de los nodos muestreados en la línea 


principal. 
Concentración de cloro (mg/l) Error absoluto promedio 
ID del nodo 
Muestreada Simulada (mg/l) Situación 
| RDTOT3126 0.86 | 114 | 0.28 | Sobreestimación | 
RDTOT3142 0.61 1.03 0.42 Sobreestimación 
| RDTOT6E420 0.76 | 0.42 | 0.34 |  Subestimación | 
DQ862 0.48 0.44 0.04 Subestimación 
DQ46 0.33 0.43 0.10 Sobreestimación 
RDTOT13 0.43 0.39 0.04 Subestimación 
CARRICAP1 0.25 0.39 0.15 Sobreestimación 
MUERTOP1 0.47 0.39 0.08 Subestimación 
QD17 0.35 0.33 0.02 Subestimación 
PEPUEBLITO 0.47 0.29 0.18 Subestimación 
QD279 0.39 0.3 0.09 Subestimación 
DQ191 0.26 0.28 0.02 Sobreestimación 
160 0.38 0.24 0.14 Subestimación 


Para obtener el resultado final de la calibración, se propone simular 
con los datos más próximos de la fecha de septiembre de 2019, que es 
cuando se tiene el modelo hidráulico calibrado. De esta manera se busca 
reducir el error por los aspectos de la variación de concentraciones ligadas 
con las características del agua desde la fuente y su efecto en la 


desinfección. Así, se propone trabajar con la concentración de inyección 
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más cercana a la fecha de septiembre de 2019, que corresponde a la 
concentración registrada por el controlador de 0.8334 mg/l (Figura 11), 


en lugar del valor promedio, manteniendo la forma del patrón de variación 
de cloro. 


Fecha (dd/mm/aaaa) 


Figura 11. Promedio diario de inyección de cloro en planta filtros. 


Además, las pruebas de frasco sugieren que la reacción del cloro con 
el agua varía de forma significativa entre septiembre de 2019 y febrero 
2020; por lo que, para obtener una representación de la información con 
mayor vinculación de las concentraciones, se usará el valor de kp obtenido 
en septiembre, que es de -1.512 días?, lo cual coincide con los muestreos 


en campo y la calibración hidráulica. 
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Después de los cambios mencionados, se realizaron tres nuevas 
pruebas de calibración. La Figura 12 compara las concentraciones 
muestreadas y simuladas en cada prueba. Las calibraciones arrojan 


concentraciones de cloro similares, lo que sugiere que se llegó a una 


calibración muy cercana a la óptima. 


ID de los nodos muestreados 
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Figura 12. Resultados del modelo en el SF812. 


Con las modificaciones realizadas respecto a la toma de datos en torno 
a septiembre de 20109, se redujo el error absoluto promedio hasta 0.080 
mg/l, y aumentó la correlación del modelo hasta 0.816, representando 
una mejoría en comparación con los resultados obtenidos durante la 


calibración inicial (Tabla 7). 


Tabla 7. Comparación de resultados de la calibración. 


Error absoluto promedio Correlación 
Núm. de 
Inicial Resultante 
prueba Inicial Resultante 
(mg/l) (mg/1) 
1 0,097 0.080 0.769 0.816 
2 0.098 0.085 0.767 0,794 
3 0.098 0.081 0.764 0.805 


La Tabla 8 presenta los coeficientes ky asignados por el CCKwDyR para 
calibrar el modelo. La Tabla 9 resume los resultados obtenidos en este 
trabajo y se compara con los logrados por otros autores. Los resultados 
dependen de la red con la que se trabajó (tamaño, nivel de 
esqueletización); la información que se tenía (información histórica o 
mediciones recientes, mediciones horarias o promedio diarias), e inclusive 


el tipo de simulación realizada (flujo permanente o no permanente). 
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Tabla 8. Coeficientes kw asignados al modelo SF812. 


13mm | ----------- IM o | o - 


13m —- RN ——--—------- 


25mm | =-1.4m/dia | -0.54 m/dia | ----------- 


a |: an | 


50mm |=-0.07 m/dia | -0.14 m/dia | ----------- 


IN | EF PE=a==<Ss> || eps | ena 


75mm IM ------- 


100mm [|=0.54 m/día | ---------=-- | ----------- 


125mm | ----------- MA o | o -- 


1i50mm |=-0.05 m/día | -0.01 m/dia  ----------- 


200mm |=0.01 m/día | -0.35 m/dia | ----------- 


250 mm A ml oo | o --- 


300 mim A | o | > 


350mm  |-=0.03 m/día | ----------- | ----------- 


600 mm  |=0.01 m/día | ---------=-- | ----------- 
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Tabla 9. Comparación de resultados. 


Error absoluto 


Autor(es) Kb (1/día) Kw (m/día) : Correlación 
medio (mg/l) 
Actual, SF812 -1.512 -0.01 a -1.50 0.080 a0.085 — 0.794 a 0.816 
Zhang, Kiene, Wable, Chan y 
“1,15 2a-2.30 | nn. oo rr O roo 
Duguet (1992) 
Vasconcelos et al. (1997) -0.23 a -1.32 -0.03 a -1.52 0.15a 0.05 0.85 a 0.98 
Elshorbagy, Abu-Qdais, 8 
sl Q AI O a 
Elsheamy (2000) 

Alcocer-Yamanaka et al. (2004) -2.3715 0.1509, -0.414 | ===. |... n---- 
Munavalli et al. (2009) -1.73 -5.24, -1.00 0.13 0.800 
Monteiro et al. (2014) -0.27 -0.04 0.0302a0.025 | —- =======---- 

Saldarriaga et al. (2014) -0.0053 A as 
Díaz-Duque (2015) -0.113 a -0.450 -0.01 a -1.29 0.120a0.107 | 0.8214 a 0.9213 

Baños-Caballero (2016) -0.104 a -0.858 -0.003 a -0.39 0.067 a 0.013 0.987 
Minaee et al. (2019) -1.73 -1,233 0.074 0.950 


El proceso de calibración descrito permite hacer una mejora del 


modelo de cloro en la red mediante un proceso de calibración incluida en 


la metodología para finalmente alcanzar un calibrado óptimo de la red. 


Uno de los aspectos importantes para lograr una aportación 


adecuada de la propuesta es la consideración del agrupamiento de 


tuberías de diámetro y rugosidad para la obtención de ky óptimo. De esta 


forma, se mantiene la diversidad de búsqueda objetiva y necesaria de la 


red, evitando un gran número de combinaciones no necesarias de valores 


de ky. Con esto se tiene una búsqueda más objetiva y enfocada a mostrar 
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resultados con bajo número de individuos discretizados que favorezcan la 
fidelidad de la calibración, como se demostró en la Tabla 2. Aplicar esta 
técnica a una red real, como la SF812, implicó no trabajar con un gran 
número de individuos debido a la gran cantidad de tuberías que tiene 
dicho sector; su esqueletización alcanzó un valor de 1 070 tuberías. Al 
agrupar las tuberías, se trabajó solo con 30 valores de kw a optimizar, 
alcanzando resultados adecuados. 


La función de aptitud propuesta es novedosa y se basa en una 
función sigmoide. Ello aporta para esta línea de investigación mediante la 
inclusión de un coeficiente de frontera para la distinción de individuos que 
obtengan un valor alto de función de aptitud respecto de aquellos que 
presenten valores bajos. En el rango de los valores bajos, su progresión 
es lenta al inicio de la técnica heurística, por lo que en las primeras 
generaciones se evita que se alcancen mínimos locales, como es el caso 
de las funciones lineales. Posteriormente, con el paso de las generaciones, 
cuando los individuos van alcanzando valores adecuados de función de 


aptitud, su progresión es más acelerada. 


De forma conjunta, la aportación de esta propuesta de investigación 
representa un proceso de calibración de modelos de calidad del agua, en 
este caso aplicado para el sector de la red SF812, que es de suma 
importancia para el organismo operador para mejorar la cloración del 
agua en Guanajuato. Con este modelo calibrado se pueden llevar a cabo 
propuestas de mejora operativa mediante los resultados de las 
simulaciones, al tener conocimiento del decaimiento del desinfectante 
dentro de la red. Por ejemplo, en la Figura 13 se muestran las 
concentraciones de cloro dentro de la red del SF812 en la hora de menor 


consumo (tiempo de residencia alto), observando que existen zonas 
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donde se pueden optimizar las concentraciones: disminuyéndola al 
comienzo de la red y aumentándola en zonas con poca recirculación de 
agua. Esto podría realizarse implementando estaciones de reinyección de 


cloro en algún punto de la red. 


Chlorine 
0.20 
0.40 
0.60 
0.80 
E malL 
E pa 
d HL 
$ 
Figura 13. Concentraciones de cloro en el SF812 en la hora de menor 
consumo. 
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Discusión 


El coeficiente kp se encuentra dentro del rango encontrado por otros 
autores. En este caso, indica un rápido decaimiento del cloro debido a una 
gran cantidad de compuestos reductores. La red se abastece de fuentes 
superficiales y se localiza en una zona minera importante (Servicio 
Geológico Mexicano, 2018), por lo que metales y otras sustancias son 
arrastrados por la cuenca hasta las presas, además de la materia orgánica 
que es común en este tipo de fuentes (Valdez, 1994). Lo anterior 
contribuye a que el valor de k» sea alto y a que no se mantenga constante 
en todos los meses del año. En el presente trabajo se empleó una 
ecuación de primer orden, por lo que en futuros proyectos se deberán 
considerar ecuaciones que consideren tanto la temperatura como el 


carbono orgánico total dentro del agua. 


Asimismo, existe una relación al incremento de los valores en el 
coeficiente kp debido a un aumento en la temperatura. Sin embargo, fue 
posible detectar valores de kp altos en muestreos en temperatura baja 
debido a la reducción de la inyección de cloro por aspectos operacionales 
del sistema de abastecimiento; se corroboró que para concentraciones 


bajas de cloro, el coeficiente de kp presenta valores elevados. 


La inyección inicial de cloro es uno de los parámetros que más 
afectan a la calibración del modelo. En la presente investigación se 
monitoreó la concentración de cloro que se inyectó a la red SF812 por 51 
días; se constató que además de la concentración constante que se realiza 


en la última etapa de la potabilización, se ve afectada por el tiempo de 
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permanencia y la temperatura variable que tiene el agua dentro del 
tanque de abastecimiento posterior a la potabilización. Se pudo obtener 
un patrón de concentración inyectado a la red, el cual fue variable y que 
finalmente, al incluirlo en la última etapa de la calibración, se mejoró en 
el modelo definitivo. Se verificó que en una mayor concentración en la 
inyección, el modelo sobreestima la concentración en los nodos más 
cercanos y que cuando se tiene una menor concentración afecta a gran 
parte del modelo, obteniendo concentraciones de cloro por debajo de los 


valores muestreados. 


Los puntos de muestreo para obtener mediciones locales del cloro 
libre residual en diversos sitios de la red se deben plantear de manera 
que sean muestreos rápidos, prácticos, que puedan realizare de manera 
fácil y sin interrumpir el suministro de agua; para ello hay que hacer una 
selección de los sitios que sean característicos de zonas de alto consumo, 
puntos alejados de la fuente y en zonas con tiempos de residencia del 
agua altos. Una buena ayuda para ello es el propio modelo hidráulico de 
la red, calibrado en periodo extendido en Epanet, para verificar tiempos 
de residencia altos y con ello tener la representatividad de la 
concentración en dichos sitios y en momentos específicos del día. En el 
trabajo de campo se alcanzó un muestro del 2.85 % dentro de la 


recomendación de la AWWA para la calibración de modelos de calidad. 


Los valores de ky suelen estar inversamente relacionado con el 
diámetro y coeficiente de rugosidad de las tuberías (Vasconcelos et al., 
1997), lo cual no ocurre con todos los ky que dan solución al modelo del 
SF812 (Tabla 9). Los coeficientes Kky propuestos por el CCKwDyR 
representan una solución matemática válida, pero no necesariamente 
corresponden a la tasa de reacción real del cloro con las paredes de las 
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tuberías e incluso a la biopelícula dentro de ésta. Como alternativa para 
mejorar aún más el modelo se podría implementar una función de aptitud 
que considere dicho aspecto y premie a los individuos que cumplan esta 


condición, y/o castigue a aquellos que no las cumplan. 


El error absoluto promedio alcanzado de 0.080 mg/l y la correlación 
de 0.816 obtenidos indican que el modelo de calidad representa el 
decaimiento de cloro de forma adecuada, pues están en el rango de 
valores obtenidos en trabajos similares anteriormente realizados. Estos 
valores podrían mejorar si se reducen los posibles errores humanos al 
momento de efectuar el muestreo, para lo cual convendría medir más de 
una vez en cada punto de monitoreo elegido. También podría ayudar un 
muestreo de datos continuo en los sitios de monitoreo, aunque vale 
recordar que es más importante la correcta distribución de los puntos 
elegidos para el muestreo, tal como lo muestran las últimas cuatro 


pruebas del análisis de sensibilidad. 


El modelo de calidad obtenido le da a SIMAPAG una buena 
representación de las concentraciones de cloro en el SF812. Implementar 
este tipo de modelos resulta una alternativa viable para cualquier 
organismo operador que busque verificar las concentraciones de cloro 


dentro y de la norma, y optimizar el uso de cloro en la red. 
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Conclusiones 


Con la implementación del CCKwDyR fue posible calibrar el modelo de 
calidad de agua del SF812. Se describió una metodología para lograr una 
correcta calibración e identificar los parámetros que lo hacen más 
sensible, a los cuales se les debe prestar mayor atención al momento de 


determinar su valor y con ello mejorar los resultados de simulación. 


El modelo calibrado representa una herramienta valiosa para los 
organismos operadores como SIMAPAG, pudiendo analizar los efectos de 
modificar la operación de la red sin necesidad de comprometer la calidad 
del agua (en los casos en los que dichas modificaciones pudieran afectar 
negativamente la calidad del agua; por ejemplo, por ajustes en la 


inyección de forma arbitraria o “ensayo y error”). 


La implementación del CCKwDyR para calibrar modelos de calidad 
de agua a través de los coeficientes Ky resulta efectiva si se determina 
con suficiente precisión el valor de parámetros necesarios: la inyección 
inicial de cloro, el coeficiente kp, el modelo hidráulico calibrado, así como 


un muestreo estratégico de las concentraciones en la red. 


Los coeficientes ky calibrados por el CCKwDyR son principalmente 
una solución matemática para calibrar el modelo de calidad. Dichos 
valores no necesariamente funcionarían para tuberías con diámetro y 
rugosidad similares, ya que en la tasa de decaimiento del cloro intervienen 
además factores como la velocidad y dirección del flujo, sustancias 


presentes en el agua, etcétera. 
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Conviene tener registro de la cantidad de cloro que se inyecta a la 
red cada hora, considerando posibles variaciones debidas a un consumo 


elevado de agua, ya que este parámetro influye mucho en la modelación. 


En caso de abastecerse de fuentes superficiales, en las cuales la 
cantidad de compuestos reductores es muy variable, conviene llevar a 
cabo pruebas de frasco constantemente para verificar que el coeficiente 


Kkp siga representando de forma correcta la tasa de decaimiento del cloro. 


Es muy importante hacer los muestreos en campo siguiendo las 
indicaciones y recomendaciones establecidas tanto por las normas 
oficiales, el organismo operador, la comisión reguladora (en México, 
Conagua) como por el fabricante del instrumento de medición utilizado, 
pues un nodo muestreado con mucha imprecisión afecta de modo grave 


a toda la calibración, principalmente a la correlación. 


El muestreo en campo, la prueba de frasco y la medición de la 
inyección de cloro se deben realizar el mismo día o al menos en días 
cercanos entre sí, evitando que las concentraciones muestreadas se 
deban a tasas de decaimiento diferentes y que esto genere errores en la 


calibración del modelo. 
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